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ABSTRACT 



The accuracy of mass estiniates by gravitational lensing using the thin- lens approx- 
imation applied to Navarro-Frenk- White mass models with a soft truncation mech- 
anism recently proposed by Baltz, Marshall and Oguri is studied. The gravitational 
lens scenario considered is the case of the inference of lens mass from the observation 
, of Einstein rings (strong lensing). It is found that the mass error incurred by the sim- 

(~| ' plifying assumption of thin lenses is below 0.5%. As a byproduct, the optimal tidal 

Qh, radius of the soft truncation mechanism is found to be at most 10 times the virial 

radius of the mass model. 
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1 INTRODUCTION 

> 

r In recent years gravitational lensing has become the main tool in use to estimate the amount and distribution of matter of 
OO astrophysical systems of various kinds and sizes, from single compact objects to clusters of galaxies. Some very wide-ranging 
predictions rest on the accuracy of such mass estimates, most notably the value of the dark energy density and, with it, the 
time since the big bang. The predictive power of gravitational lensing as a tool to study individual deep objects also hinges 
. heavily on the accuracy of the method. Consequently, any studies of the degree of accuracy of the method are naturally 
' warranted. In this paper we describe a study that provides a quantitative measure of the degree of accuracy of the method 
^ : in a particular context. Together with those of a previous study of ours ( |Kling fc Frittelli 2008 [ ), these results are among the 
• • first quantitative estimates of the accuracy of the thin-lens approximation applied to gravitational lensing systems. 
^ ^ , Our conceptualization of a gravitational lensing system is as follows: There is a physical "lens" that gravitationally 

' disturbs the path of lightrays from a source to the observer. The theoretical framework for this situation would be a spacetime 
$^ : with a metric that solves the Einstein equations with a mass density source that describes the lens. The lightrays from the 
source to the observer would be null geodesies connecting both points. Even though the gravitational lensing event in itself 
(the null rays) is very amenable to a straightforward computational strategy, there is an intrinsic difficulty with this theoretical 
framework: the metric. This is determined by the mass density of the lens, also known as the density profile, and this, at this 
time, is still subject to modeling. 

The conventional approach consists in making use of the thin-lens approximation to calculate the bending of lightrays. 
The conventional approach also makes use of the density profile assumed for the physical system. However, the quantities 
of interest in the conventional approach require a set of manipulations on the density profile that differ significantly from 
those needed in the theoretical framework. One significant difference resides in the fact that most quantities of interest in 
the conventional approach require the mass density integrated along the line of sight, a quantity referred to as the projected 
mass. 

The conditions that a mass density must meet as a function in space in order to be useful for the conventional approach 
are much weaker than what is needed for the theoretical framework. The conventional approach admits density profiles that 
are not integrable over all space, namely, that have slow fall-off rates (slower than r~^). Admittedly, such density profiles 
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are unphysical in the sense that they do not represent the isolated systems that they intend to model. Nevertheless, they 
continue to be used because they are known to be a close approximation to the underlying physical system within a finite 
region centered on the lens. The question arises as to how to "terminate" such unphysical density profiles in the region where 
the approximation is not valid anymore (i.e., far from the lens). 

The simplest way to terminate an unphysical density profile is by "hard" truncation, that is: by setting the profile 
discontinuously to zero beyond some radius. Arguably, this mechanism results in a density profile with questionable physical 
sense, especially because the truncation radius is largely arbitrary but also because discontinuous truncation introduces 
significant spurious magnification of images. However, the underlying system may be better approximated by such a drastically 
truncated density profile than by a non-integrable one, so the method is not devoid of merit. 

Alternatively to hard truncations, various "soft" truncation mechanisms have been proposed, which normally consist of 
multiplicative factors that speed up the fall-off rate of the density profile in order to ensure a finite total mass. 

In the present work we investigate how well the conventional approach to gravitational lensing estimates the mass of a 
lens, compared to the theoretical prediction. There are several subtleties associated with this seemingly straightforward task. 
To start with, the physical system itself (the actual density profile of a cluster in all space) is not known. What is known is 
a good representation of the density profile within a certain radius, referred to as the virial radius. Since the "true" density 
profile of the lens must be integrable over all space, we take the view that a profile with a finite total mass is the "true" 
system, and any other profile with a fall-off rate slower than r^^ is an approximation to the true system, whose degree of 
accuracy must be quantified. In this view, a profile with an infinite total mass takes second seat to its truncated version, 
regardless of the truncation mechanism. This may sound contrary to prevailing practice, but is absolutely necessary to our 
goal because, in the theoretical framework, non-physical profiles result in inaccurate lightrays. 

Taking a given truncated profile as the "true" physical system, there are thus two different approximations to quantify 
in relation to the "true" mass of the system as predicted by the theoretical approach: the error introduced by the thin-lens 
approximation and the further error introduced by the use of a profile that is integrable along the line of sight but which 
results in an infinite total mass. In this paper we refer to the latter as the error introduced by "removal of truncation", for 
lack of a better phrase. Both errors are quantified in the present work within a given family of density profiles. 

One can, and probably should, think of the truncation mechanism itself as an approximation: soft versus hard, for 
instance, or different kinds of soft truncation, all applied to the same model. In that case, it makes sense to quantify the error 
introduced by the truncation mechanism as a third source of inaccuracy in the mass estimates. We are able to provide an 
estimate of the truncation mechanism error by comparing mass predictions using soft and hard truncation schemes for the 
same class of density profiles. 

The gravitational lens scenario considered is the case of the inference of lens mass from the observation of Einstein rings 
(strong lensing) . The family of density profiles of interest is described in Section [J] The mass prediction by the conventional 
approach is discussed in Section |3] The calculation of the "true" mass as predicted by the theoretical framework is found in 
Section [l] The results of our study are found in Section [5] Section |6l We conclude in Section [7] with some remarks of interest 
that follow from our study. 



2 DARK MATTER HALOES 

Numerical modeling of dark matter haloes predicts a density profile of the form ([Navarro et al. 1997|l 

p{r) _ 5^ 



(1) 



hereafter referred to as the NFW model, where pcrit is the critical density (given by 3H{z)^ /SivG in terms of the gravitational 
constant and Bubble's constant at the lens redshift), is a scale radius defined as the peak of r^p(r), and Sc is a characteristic 
density contrast. A "virial" radius is defined as the radius r2oo of the sphere with mean density equal to 200pcrit, and is used 
to identify collapsed particle systems in numerical simulations. The ratio r2oo/rs = c is referred to as the halo's concentration. 
By definition of the virial radius one has, as a consequence, 

^ 200 

3 log(l + c) - c/(l + c)- ^ ' 

Falling off as r^^, the NFW profile is not integrable over all space, which is interpreted as an indication that the model 
fails at large radius. Hence, a truncation mechanism is needed for long-range applications, such as gravitational lensing. Two 
truncation mechanisms for the NFW model have been used: hard truncation (|Takada fc Jain 20031 |Kling fc Frittelli 2008[ ) 
and smooth truncation (|Baltz et al. 2009|) . The hard-truncation mechanism consists of discontinuously terminating the model 
by assuming that the mass density is identically vanishing outside a given radius. In the lack of good knowledge pertaining to 
halo edges, the choice of termination radius is largely arbitrary. The common choice has been to terminate the NFW model 
at the virial radius. 
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The smooth truncation mechanism of Baltz at al. (|2009[) consists of multiplying the NFW profile by a factor that falls 
off at least like r~^, being close to 1 within small values of r. This mechanism yields a sufficiently fast decay to ensure a finite 
total mass (at least r~^), while preserving the original profile within a finite radius, within some accuracy. The smoothly 
truncated profile of Baltz et al. I|2009[l is 

pAn = 77 — (3) 

where n is a positive integer and rt is a parameter referred to as the tidal radius. For the purposes of the current investigation, 
we restrict attention to the case n = f, hereafter referred to as the truncated profile. For notational convenience we use po 
for the NFW model, since clearly the expression ([Sjl reduces to the NFW model ^ for n = 0. If we use a rescaled coordinate 
X = r/vs, the virial radius would be located at x = c and the tidal radius would be at a; = r = rt/rs- At any given radius x, 
the truncated profile pi differs from the NFW profile po by 

Pijx) - po{x) ^ x^ 

po{x) x^ + T^- ^ ' 

As a consequence, one way to interpret the tidal radius is as the radius at which the profiles differ by 50%. At the virial radius 
the profiles differ by 

Pi(c)-po(c) 1 



(5) 



Large values of the tidal radius r relative to the virial radius thus ensure that the truncated profile remains sufficiently close 
to the NFW profile within the virial radius. For the profiles to differ by 10% or less at any point within the virial radius, the 
tidal radius needs to be chosen as r = 3c or larger. For this range of values of r, the virial mass (the mass contained within 
the virial radius) of the truncated profile differs from the virial mass of the NFW profile by less than 3%. 

On the other hand, for large values of the tidal radius the truncated profile approaches the NFW profile, acquiring the 
undesirable features of the NFW model, such as a very large unrealistic total mass. The optimal value of r would thus be 
the smallest possible value that yields the desired accuracy within the virial radius, which is somewhat arbitrary itself. Baltz 
et al. (2009 ) justify the choice of r = 2c on the basis that for this tidal radius the virial mass of the truncated profile differs 
from the NFW virial mass by only 6%. With this choice, however, the profiles themselves differ by up to 20% at the virial 
radius, a difference that may be considered significant. At any rate, the tidal radius should not be viewed as an independent 
parameter but should be considered a function of the concentration parameter. With this assumption, the models have two 
free parameters: the concentration parameter c and the scale radius Vs- 



3 MASS ESTIMATES BY THIN-LENS GRAVITATIONAL LENSING 

We are interested in estimating the error in the prediction of the mass of a cluster by the observation of Einstein rings 
produced by the perfect alignment of a source of light behind a massive spherically symmetric deflector along the line of sight. 
In the usual manner (|Schneider et al. 1992|) . the radius s of the Einstein ring on the deflector's plane (perpendicular to the 
line of sight) at a distance Di from us is the root of the lens equation: 

where M„ (s) is the projected mass within a distance s to the center of the lens: 

M^{s) = 27V s'ds' J p„ (^^/s'^ + p'j M. (7) 
The constant Ecr is the characteristic surface mass density, given by 

where Ds is the distance to the source, Dis is the distance between the lens and the source, G is the universal gravitational 
constant, and c is the speed of light, not to be confused with the NFW concentration parameter. 

For a given conflguration of source and lens, different spherically symmetric models p,i will predict different sizes of 
Einstein rings by virtue of the difference in their projected masses Mn{s). In principle, given a known Einstein ring and 
known source and lens locations, equation ((6| determines the concentration co or ci of the NFW model or the truncated 
model, respectively, as functions of the scale radius . Different concentrations co and ci would lead to different predicted 
virial masses. In practice, we flx the value of the scale radius, the Einstein ring's angle and the distances to the lens and 
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the source, and determine the concentration by numerically solving equation ((6| for c„ to high accuracy using a standard 
bisection algorithm. 

The value of the concentration parameter determines the virial mass of the model M^o^ by; 



Ml°° = 47r / r''p„{r)dT (9) 







which, in the case of n = 0, is equivalent to Afo°° = CoTr^pcrit?"? 



4 MASS ESTIMATES BY NUMERICAL INTEGRATION OF THE LIGHTRAYS 

The standard gravitational lensing formalism, leading to Eq.([6|, assumes that the lightrays from the source to the observer 
are straight lines with sharp bending at the plane containing the deflector. The bending and all its consequences are obtained 
by "projecting" the mass density along the line of sight, namely, by integrating the mass density along the line of sight. 
This formalism is shadowed by a general relativistic formalism in which the lightrays from the source to the observer curve 
according to the geodesic equation of the metric of spacetime. The spacetime in this case is a lens with a weak gravitational 
field located on a expanding cosmological background. The lens is described by the truncated model pi. 

There is some leeway as to how to put together the metrics of an isolated lens and that of a Friedman-Robertson- Walker 
cosmology, in the sense that, to our knowledge, a rigorous and complete derivation of the Einstein equations valid for the 
particular case of an isolated lens on a cosmological background is still lacking. The Einstein equations matter because the 
motion of the spacetime affects the matter density, which affects the metric in turn. So long as the rigorous results remain 
outstanding, plausible (if unproven) descriptions of the behavior of the matter in the lens may need to be used. A widely 
accepted notion, supported by preliminary results, is that the matter distribution in a cluster of galaxies does not partake 
of the general expansion (namely, clusters do not expand internally). This can be modeled with a metric that approaches 
the standard cosmology at large distances from the lens, but which, in the vicinity of the lens itself, is indistinguishable from 
an isolated lens in flat space. Such a metric would be of the form 

ds^ = {l + 2^p{r,t))dt^ - a{tf{l - 2ip{r,t)){dr^ + r^(de^ + sin^ (10) 

where if represents the non-expanding isolated lens centered at r = 0. Wherever the function ip is not negligibly small, the 
space part of the metric would be equivalent to that of an isolated lens on flat space if ip{r, t) is the Newtonian potential of 
the lens in proper coordinates, namely ip{r,t) — U{a{t)r) with \7^U — AiiGpi. By requiring the solution of this equation to 
vanish at infinity and be finite at the origin, the Newtonian potential U of the mass distribution pi is found to be 



Gm 



(l+r2 



arctan(./.) ( - - . - 2-) + In ) ' + " ^ In . 



(11) 



with X = r /ts. Ifere m = 47rr3 5cPcrit. 

The function tf{r,t) is obtained by the substitution r — > a{t)r in Eq. and is, in principle, a time-varying function. 

However, the time variation is extremely slow in the conditions of observable lensing, where the source and observer are very 
far away from the lens. Since ip is relevant only in the vicinity of the lens, one could substitute the time-varying 
a{t) with a constant a{ti) with ti being a representative time at which the lightray being considered reaches 
the vicinity of the lens. A posteriori, we find that the results of the numerical calculations that follow are not significantly 
affected by the choice of flxing the value of a{t) to a constant, so we follow common practice and set 

p>{r,t)^U{a{U)r) (12) 

for use in the metric pO|l . This has the additional advantage that in these conditions the metric (|10|l is conformally static, 
as assumed in Schneider et al (|1992p . and thus the well known phenomena associated with lenses follow as per that text, 
including Eq. 6. 

In the standard cosmology, the expansion parameter a{t) is given by 

with Hq = 70 km/s/Mpc, f7m = 0.3 and f^A = 0.7 = 1 - f7m. 

Since the metric is spherically symmetric, we can restrict attention to the equatorial plane without loss of generality. The 
null geodesic equations are equivalent to the the Euler-Lagrange equations of the Lagrangian 

C=(l + 2ip)P - a^(t)(l - 2p){r'' + r^^} = 0, (14) 



© 0000 RAS, MNRAS 000, 000-000 



Thin lens accuracy 5 



The Lagrangian is set to zero because the geodesies are null. The Euler-Lagrange equations are equivalent to five first-order 
ODEs, which we can write as 



Vt 



davt „ .d(p 

Vt = --r- 2(1 - 2<p)—VtVr 

at a or 

2vtVr d-CL _^ ^<P d(p 2 

^ a dt dr ^ 

26^ dif,^ „ , 



4> = -^(1 + 2^). (15) 

These equations are accurate to first order in tp. The parameter & is a constant of integration and is related to the "observation 
angle" at the observer, or the angle between the lightray and the optical (radial) axis connecting the observer to the lens as 
explained in Kling & FritteUi IpOOS)) : 

sin 6*063 = Sobs = —(1 + 2v3o), (16) 
ro 

where the potential ipQ, is evaluated at the observer position, ro. 

By integrating the equations of motion from a given point, for a specified set of parameters (r^, c) describing the matter 
density, one finds a specific light ray passing through that point. Our scheme for a mass prediction is based on lightrays 
making Einstein rings, for which the source, lens and observer are coUinear and are at known distances from each other. So 
the beginning and end point of the null geodesies being sought are known. In these conditions, for every observed Einstein 
ring Oobs, as per Eg. 1161 there is an associated set of values of (rs,c) that are consistent with the observed angle as well as the 
fixed configuration of lens, source and observer. As in the case of the prediction by the thin-lens gravitational lens formalism, 
in practice we fix Ts and find c using Newton's method and a ray shooting technique. The ray shooting technique consists 
essentially of integrating rays for a given value of c and seeing where they land with respect to the end point aimed at, either 
in front of it or behind it, subsequently adjusting c accordingly, until the ray integrated lands sufficiently close to the target, 
within a specified accuracy. In practice, the integration is done from observer to source. 

For the purposes of integrating the null geodesies, we re-scale the time and radial coordinates by the age of the universe. 
The equations of motion are integrated using an adaptive step-size Runge-Kutta-Fehlberg 4-5 method based on the imple- 
mentation in (|Press et al. 1995|l . The numerical error is at least 500 times smaller than the difference between the result by 
null geodesic integration and the prediction by the thin-lens formalism. 

The value of c found by this method is considered to be the "true" value of the model and is used to predict the "true" 
virial mass of the model M^"" by Eq. ([9]) with n = 1, namely: 

M^"'' = 47r / rVi(r)dr. (17) 



5 COMPARISON OF THE THIN-LENS MASS PREDICTION RELATIVE TO THE NULL-GEODESIC 
PREDICTION 

For our first study we arbitrarily choose a series of five reasonable lensing systems with varying redshifts, and calculate 
the concentration parameters (c, ci,co) by the method of null-geodesies, thin-lens with truncation and thin-lens with no 
truncation, respectively. We do this for r = 2c and for r = 3c separately. In all cases, the scale radius is arbitrarily assumed 
to be rs = 0.25 Mpc. The results are shown in Tables [T] and [2] 

In all cases, the concentration parameter ci is slightly higher than the true value c, whereas cq is closer to and generally 
lower than c. However, the reader should keep in mind that when the virial mass is used as the parameter for 
comparison, the relative values of cq and ci are not significant in themselves. The tables show for each case the 
errors in the virial masses of the truncated and non-truncated thin lens models relative to the true model (which is truncated 
with no assumption of thin lenses). 

We can verify that, whether truncation is applied or not, the use of the thin-lens approximation overestimates the virial 
mass of the system, as the relative errors are negative in all cases. 

The sizes of the two separate simplifying assumptions (thin lens and removal of truncation) are revealed to be very 
different. The error incurred by the sole simplifying assumption of thin lenses, displayed in the column headed by 1 — 
Mf°/M^°°, is generally between 0.3% and 0.7% for r = 2c, and between 0.2% and 0.4% for r = 3c. 
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By contrast, the error incurred by further removing the truncation in order to simplify the estimation of the mass is about 
6% for r = 2c and 3% for r = 3c. This is the dominant source of error in the virial mass. Keeping in mind that "truncation" 
in this model refers merely to a fast decay, this is not surprising, as different decay rates lead to significantly different profiles 
within the virial radius, as shown in Figure [1] 

For our second study, we consider three observed systems with symmetrical arcs where one can use an Einstein ring 
analysis to determine parameters of the system with confidence. The systems are RXJ1347-1145, the brightest source in the 
ROSAT all sky survey (|Sahu et al. 1998|) : MS2137-23, a relatively structure- free, four-arc system l|Gavazzi et al. 2003|) : and 
cluster A of the high-redshift cluster RCS2319-I-00 l|Gilbank et al. 2008|) . Tabled gives the redshift and arc distribution details 
of these clusters. 

For each one of these three systems, the Einstein rings are known. Modeling each system as a truncated pi density profile 
with r = 3c leads to two unknown parameters for each system: the value of the concentration parameter c and the value of 
the scale radius Ts. The left panel of Figure [2] shows the contour plot in the — c plane for the Einstein rings of these three 
systems. For this plot, was varied arbitrarily and c is the true concentration parameter determined using the null-geodesic 
integration formalism. In the right panel of Fig. (2] we show the error in the virial mass 

l-Mr/M^"". As the errors are 

negative, we verify that the thin-lens approximation overestimates the mass of each system. For reasonable values of Ts and 
c, the errors are of the order of 0.5%. 

Figure [3] shows the corresponding results for the case that the modeling is done with r = 2c. One can verify that even 
though the concentration parameters do not differ significantly from the ones in the case of r = 3c, the accuracy of the mass 
prediction is reduced by a factor of 2, as the maximum relative error in the mass is of the order of 10%. 

6 COMPARISON OF TRUNCATION MECHANISMS 

As indicated in the introduction, one of the fundamental problems of gravitational lensing is the level of uncertainty in the 
underlying physical model of the lens irrespective of the approximation method used to make predictions. Here we address 
some of the consequences of this uncertainty by calculating the difference between the virial masses predicted exactly (with 
null geodesies) for two different models assumed to represent the same physical lens. 

The two models differ in the truncation mechanism applied to an NFW lens: soft versus hard truncation. We use a density 
profile pi with t = 3c as a softly truncated NFW lens with rs = 0.25 Mpc. The hard truncation model is a density profile 
po with Ts = 0.25 Mpc within the virial radius, and zero outside of the virial radius. The hard truncation model is discussed 
extensively in our previous work ( |Kling &i Frittelli 2008[ ), where the gravitational potential needed for the integration of the 
null geodesies is found (Equation (34) in ( |Kling fc Frittelh 20 08')). 

Table[4]shows the result of this study. The predicted virial masses for both models differ within 4%. This is an independent 
indication that the truncation mechanism itself is much more important than the simplifying assumption of thin lenses, and is 
the dominant source of error in mass estimates by gravitational lensing, by an order of magnitude. Perhaps counter-intuitively, 
the difference in profiles outside of the virial radius (as seen in Figure [l| is a significant source of inaccuracy in the prediction 
of virial mass by actual light paths traveled. 

One puzzling result is that the hard truncation model consistently overestimates the mass with respect to the soft 
truncation model (as indicated by the negative values in the last column in Table [4]). With both profiles being almost identical 
within the virial radius, one would expect that the soft truncation mechanism (just be virtue of having more total mass) 
would produce larger bending and would lead to a higher estimate for the virial mass. 

7 DISCUSSION 

One of the main results obtained through our present study is that given a reasonable physical mass profile, such as pi, the 
use of the thin-lens approximation as a simplifying assumption to predict the virial mass of the system leads to very small 
errors, which we found to be of order 0.5%. This can be interpreted as a validation of the use of the thin-lens approximation 
in the present context and, by extension, as an indication of the validity of the thin-lens approximation in generic physically 
reasonable systems. 

A second result found is that, for this family of models, given a reasonable physical mass profile, such as pi, removing 
the truncation as a further simplifying assumption (such as using the NFW) in addition to the thin-lens approximation leads 
to a significant loss of accuracy in the prediction of the virial mass of the system in the cases r = 2c and r = 3c. The error 
introduced by removing the truncation is larger by a factor of at least 10 over the error introduced by the use of the thin-lens 
approximation. 

For practitioners of the thin-lens approximation, one of the main concerns with the soft truncation mechanism studied 
here is the built-in r-dependent discrepancy in the mass profile. In effect, systems modeled with profiles pi with different 
values of r/c can vary significantly. Because the truncation is effected by means of a modification of the decay rate, the profiles 
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diverge over a significant range of distances. What is important, liowever, is tfiat tfie r— profiles are reasonably close to the 
NFW profile po within the virial radius. This is because the NFW profile is known to be accurate up to the virial radius, so 
any loss of accuracy within the virial radius is a disadvantage of the r— dependent model. Baltz et al H2009[) propose that a 
discrepancy of about 20% in the mass profile between pi and po is acceptable, as the penalty for working with a small tidal 
parameter r — 2c. 

However, larger values of the tidal parameter lead to better accuracy in the mass profile within the virial radius. We have 
argued and demonstrated that a slightly higher value r = 3c cuts the error in both the virial mass and density profile by half. 
Raising the value of t relative to c increases the accuracy of the profile, in principle arbitrarily. However, a profile pi for high 
r can approach po too accurately to be useful, if the predicted virial mass M^"" approaches Mq"" more accurately than the 
limiting accuracy of the thin-lens approximation itself. In other words, since the thin-lens approximation itself carries a base 
accuracy of 0.5% in the virial mass prediction, there is no use in fine-tuning the density profile to lead to an accuracy better 
than 0.5% in the virial mass. 

Our present study thus allows us to formulate a criterion for an optimal value of the tidal parameter r relative to c: r/c 
should be as large as possible but larger than required to achieve an accuracy of about 0.5% in the virial mass Mi°° relative 
to the NFW virial mass Mq"". We can come up with a rough estimation by noticing that 20% accuracy in the profile leads to 
6% accuracy in the virial mass, whereas 10% accuracy in the profile leads to 3% accuracy in the virial mass. The gross trend 
would indicate that 0.1% accuracy in the profile would achieve roughly 0.3% accuracy in the virial mass. By equation (5] a 
relative error of no less than 1% in the profile leads to a tidal radius no larger than r ~ 10c. A value of the order of r = 10c 
thus would guarantee that the NFW profile is modeled with pi within the virial radius to the limit of accuracy of the thin 
lens approximation, while ensuring that the model has a finite total mass. 

This rough estimate is validated by Table where the errors in the virial masses Mj and Mq"" relative to the true 
virial mass M^"*^ are shown in the last two columns. The column headed by 1 — Mf^/M^"" represents the error incurred 
by the use of the thin lens approximation, and is generally of order of 0.2%. The last column, headed by 1 — M^"'^ / M^"'^ , 
represents the combined error of the thin lens approximation and lack of truncation, and is generally of order 0.4%. The 
difference between the values in both columns is the error introduced by the truncation mechanism alone, and is of order 
0.2%, namely, comparable to the thin-lens approximation error, as anticipated. We claim that for thin-lens practitioners, using 
T much greater than 10c would be inefficient, as it would lead to a truncation error smaller than the error inherently built 
into the thin-lens approximation itself. 

Further validation of this claim is provided by Figure |4l where the relative virial mass errors 1 - MfO/M^oo and 
are calculated for values of r = jc with j — 1,2, ...15, for the three observed systems of Table [5] One can 
see that in all three cases the thin-lens errors hover around 0.2%, whereas the combined truncation/thin- lens error drops 
dramatically from roughly 2% at r = 3c to less than 0.5% at r = 8c, reaching within a factor of two of the thin-lens error at 
T = lOc, indicating that the truncation error alone is of the same size as the thin-lens error at r = 10c, as anticipated. 

Other features to note from inspection of Tables [T] [2] and [3] are as follows. The larger the predicted mass, the more accurate 
the thin-lens prediction, with or without truncation. Nevertheless, the non-truncated prediction improves only marginally in 
accuracy, in contrast to the truncated prediction, which improves significantly. 
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T = 3c 


System 




r200 (Mpc) 


c 




Co 


M200 (IOI^Mq) 


1-Mf0/M200 


1 - M^oo/M^oo 




10.0 


2.05 


8.201 


8.210 


8.204 


1.14 


-0.0033 


-0.0289 


zi = 0.2 


17.5 


2.29 


9.153 


9.161 


9.154 


1.60 


-0.0028 


-0.0274 




25.0 


2.48 


9.939 


9.948 


9.940 


2.04 


u.uu^u 


-0.0264 


Zs = 0.4 


32.5 


2.66 


10.639 


10.648 


10.639 


2.50 


-0.0025 


-0.0257 




40.0 


2.82 


11.282 


11.291 


11.281 


2.98 


-0.0025 


-0.0251 




10.0 


1.69 


6.756 


6.763 


6.755 


0.72 


-0.0031 


-0.0289 


zi = 0.3 


17.5 


1.91 


7.658 


7.664 


7.655 


1.04 


-0.0025 


-0.0270 




25.0 


2.10 


8.415 


8.421 


8.410 


1.38 


-fi nn99 


-0.0257 


Zs = 0.8 


32.5 


2.27 


9.093 


9.100 


9.087 


1.74 


-0.0021 


-0.0248 




40.0 


2.43 


9.720 


9.727 


9.712 


2.12 


-0.0020 


-0.0240 




5.0 


1.34 


5.376 


5.384 


5.377 


0.36 


-0.0045 


-0.0319 


zi = 0.3 


9.0 


1.49 


5.972 


5.979 


5.971 


0.50 


-0.0035 


-0.0297 




1 n 


1 fii 










n nn9Q 

-u.uuzy 


-U.UZOO 


Zs = 1.5 


17.0 


1.72 


6.874 


6.880 


6.870 


0.76 


-0.0026 


-0.0273 




21.0 


1.81 


7.258 


7.264 


7.253 


0.88 


-0.0024 


-0.0265 




10.0 


2.05 


8.198 


8.207 


8.199 


1.82 


-0.0032 


-0.0280 


zi = 0.6 


17.5 


2.37 


9.466 


9.475 


9.464 


2.80 


-0.0029 


-0.0262 




25.0 


2.64 


10.546 


10.556 


10.544 


3.86 


-0 0029 


-0.0251 


Zs = 0.8 


32.5 


2.88 


11.521 


11.533 


11.518 


5.04 


-0.0030 


-0.0243 




40.0 


3.11 


12.426 


12.439 


12.422 


6.32 


-0.0032 


-0.0237 




5.0 


1.30 


5.205 


5.214 


5.205 


0.74 


-0.0049 


-0.0313 


zi = 1.0 


9.0 


1.48 


5.933 


5.941 


5.930 


1.10 


-0.0038 


-0.0286 




13.0 


1.63 


6.538 


6.546 


6.533 


1.48 


-0.0033 


-0.0269 


Zs = 1.5 


17.0 


1.77 


7.078 


7.086 


7.071 


1.88 


-0.0030 


-0.0257 




21.0 


1.89 


7.576 


7.584 


7.567 


2.30 


-0.0029 


-0.0247 



Table 1. The concentration parameters and relative errors in the virial mass for a set of hypothetical lensing systems assuming that the 
true system has a density profile pi with a scale radius of rs = 0.25 Mpc and t = 3c. We give the true c value, and values obtained by 
thin lens approximations using a density profile pi (ci) and an NFW profile po (co)- The virial radius quoted, r2oo> is the one calculated 
using the true c value. The true mass, M^"" is given in units of 10^^ solar masses. 
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T = 2c 


System 




r200 (Mpc) 


c 


Cl 


Co 


M200 (IOI^Mq) 


1-Mf0/M200 


1 - M^oo/M^oo 




10.0 


2.05 


8.202 


8.216 


8.204 


1.12 


-0.0049 


-0.0611 


zi = 0.2 


17.5 


2.29 


9.156 


9.168 


9.154 


1.54 


-0.0040 


-0.0579 




25.0 


2.49 


9.944 


9.956 


9.940 


1.98 


u.uuoo 


u.uuuo 


Zs = 0.4 


32.5 


2.66 


10.647 


10.656 


10.639 


2.44 


-0.0033 


-0.0541 




40.0 


2.82 


11.289 


11.300 


11.281 


2.90 


-0.0031 


-0.0528 




10.0 


1.69 


6.757 


6.770 


6.755 


0.70 


-0.0054 


-0.0626 


zi = 0.3 


17.5 


1.92 


7.662 


7.673 


7.655 


1.00 


-0.0041 


-0.0583 




25.0 


2.11 


8.421 


8.431 


8.410 


1.34 


-fl 00*?^ 




Zs = 0.8 


32.5 


2.28 


9.102 


9.111 


9.087 


1.70 


-0.0031 


-0.0533 




40.0 


2.43 


9.731 


9.740 


9.712 


2.08 


-0.0029 


-0.0515 




5.0 


1.34 


5.374 


5.389 


5.377 


0.34 


-0.0087 


-0.0696 


zi = 0.3 


9.0 


1.49 


5.973 


5.986 


5.971 


0.48 


-0.0066 


-0.0649 




1 n 


1 fii 








u.uu 




-u.uu±o 


Zs = 1.5 


17.0 


1.72 


6.878 


6.889 


6.870 


0.72 


-0.0047 


-0.0595 




21.0 


1.82 


7.264 


7.274 


7.253 


0.86 


-0.0042 


-0.0577 




10.0 


2.05 


8.202 


8.214 


8.199 


1.76 


-0.0047 


-0.0591 


Zi = 0.6 


17.5 


2.37 


9.472 


9.485 


9.464 


2.72 


-0.0039 


-0.0551 




25.0 


2.64 


10.556 


10.568 


10.544 


3.76 


-0 0036 


-0 0525 


Zs = 0.8 


32.5 


2.88 


11.533 


11.547 


11.518 


4.92 


-0.0035 


-0.0505 




40.0 


3.11 


12.440 


12.455 


12.422 


6.18 


-0.0036 


-0.0489 




5.0 


1.30 


5.206 


5.221 


5.205 


0.72 


-0.0091 


-0.0680 


zi = 1.0 


9.0 


1.48 


5.937 


5.950 


5.930 


1.06 


-0.0067 


-0.0622 




13.0 


1.64 


6.545 


6.557 


6.533 


1.44 


-0.0055 


-0.0584 


Za = 1.5 


17.0 


1.77 


7.087 


7.099 


7.071 


1.82 


-0.0048 


-0.0555 




21.0 


1.90 


7.587 


7.598 


7.567 


2.24 


-0.0044 


-0.0532 



Table 2. The concentration parameters and relative errors in the virial mass for a set of hypothetical lensing systems assuming that the 
true system has a density profile pi with a scale radius of rs = 0.25 Mpc and t = 2c. We give the true c value, and values obtained by 
thin lens approximations using a density profile pi (ci) and an NFW profile po (co). The virial radius quoted, r2oo> is the one calculated 
using the true c value. The true mass, M^"" is given in units of W^^ solar masses. 
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T = lOc 


System 




r200 (Mpc) 


c 


Cl 


Co 


M200 (IOI^Mq) 


1-Mf0/M200 


1 - M^oo/M^oo 




10.0 


2.05 


8.199 


8.205 


8.204 


1.18 


-0.0022 


-0.0045 


zi = 0.2 


17.5 


2.29 


9.149 


9.155 


9.154 


1.64 


-0.0020 


-0.0042 




25.0 


2.48 


9.934 


9.941 


9.940 


2.08 


-fi nn9n 

u.uu^u 


-0.0041 


Zs = 0.4 


32.5 


2.66 


10.633 


10.640 


10.639 


2.56 


-0.0020 


-0.0041 




40.0 


2.82 


11.274 


11.282 


11.281 


3.06 


-0.0021 


-0.0041 




10.0 


1.69 


6.753 


6.756 


6.755 


0.74 


-0.0014 


-0.0037 


zi = 0.3 


17.5 


1.91 


7.652 


7.656 


7.655 


1.06 


-0.0013 


-0.0034 




25.0 


2.10 


8.408 


8.411 


8.410 


1.42 


-0.0013 


-0.0034 


Zs = 0.8 


32.5 


2.27 


9.084 


9.089 


9.087 


1.78 


-0.0014 


-0.0033 




40.0 


2.43 


9.710 


9.714 


9.712 


2.08 


-0.0014 


-0.0033 




5.0 


1.34 


5.375 


5.378 


5.377 


0.36 


-0.0015 


-0.0038 


zi = 0.3 


9.0 


1.49 


5.970 


5.972 


5.971 


0.50 


-0.0012 


-0.0035 




1 n 


1 fii 


6 449 








n nni i 

-U.UUl -L 


-u.uuoo 


Zs = 1.5 


17.0 


1.72 


6.869 


6.871 


6.870 


0.76 


-0.0011 


-0.0032 




21.0 


1.81 


7.252 


7.254 


7.253 


0.90 


-0.0011 


-0.0031 




10.0 


2.05 


8.193 


8.200 


8.199 


1.86 


-0.0023 


-0.0045 


zi = 0.6 


17.5 


2.36 


9.458 


9.466 


9.464 


2.86 


-0.0023 


-0.0044 




25.0 


2.63 


10.536 


10.545 


10.543 


3.94 


-0 0025 


-0.0044 


Zs = 0.8 


32.5 


2.88 


11.510 


11.520 


11.518 


5.14 


-0.0027 


-0.0046 




40.0 


3.10 


12.412 


12.424 


12.422 


6.44 


-0.0030 


-0.0047 




5.0 


1.30 


5.203 


5.206 


5.205 


0.76 


-0.0020 


-0.0042 


zi = 1.0 


9.0 


1.48 


5.928 


5.932 


5.930 


1.12 


-0.0018 


-0.0039 




13.0 


1.63 


6.531 


6.535 


6.533 


1.50 


-0.0018 


-0.0037 


Zs = 1.5 


17.0 


1.77 


7.069 


7.073 


7.071 


1.92 


-0.0018 


-0.0037 




21.0 


1.89 


7.565 


7.569 


7.567 


2.34 


-0.0019 


-0.0037 



Table 3. The concentration parameters and relative errors in the virial mass for a set of hypothetical lensing systems assuming that the 
true system has a density profile pi with a scale radius of rs = 0.25 Mpc and t = 10c. We give the true c value, and values obtained by 
thin lens approximations using a density profile pi (ci) and an NFW profile po (co). The virial radius quoted, r2oo> is the noe calculated 
using the true c value. The true mass, M^"" is given in units of 10^^ solar masses. 
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rs = 0.25 Mpc 



System 




r|°ft (Mpc) 


Csoft 


Chard 


M200 (lO^Mo) 


1 _ m200^/M200 
hard' sort 


zi = 0.2 
Zs = 0.4 


10.0 
17.5 
25.0 
32.5 
40.0 


2.05 
2.29 
2.48 
2.66 
2.82 


8.201 
9.153 
9.939 
10.639 
11.282 


8.209 
9.161 
9.947 
10.647 
11.290 


1.15 
1.59 
2.04 
2.51 
2.99 


-0.0308 
-0.0296 
-0.0287 
-0.0281 
-0.0276 


zi = 0.3 
Zs = 0.8 


10.0 
17.5 
25.0 

32.5 
40.0 


1.69 
1.91 
2.10 

2.27 
2.43 


6.756 
7.658 
8.415 
9.093 
9.720 


6.766 
7.668 
8.425 
9.104 
9.731 


0.71 
1.04 
1.38 
1.74 
2.13 


-0.0340 
-0.0324 
-0.0314 

-0.0306 
-0.0300 


zi = 0.3 
Zs = 1.5 


5.0 
9.0 
13.0 
17.0 
21.0 


1.34 

1.49 
1.61 
1.72 
1.81 


5.376 

5.972 
6.453 
6.874 
7.258 


5.388 

5.984 
6.464 
6.886 
7.270 


0.36 
0.49 
0.62 
0.75 
0.88 


-0.0380 

-0.0363 
-0.0351 
-0.0343 
-0.0336 


zi = 0.6 
Zs = 0.8 


10.0 
17.5 
25.0 
32.5 
40.0 


2.05 
2.37 
2.64 
2.88 
3.11 


8.198 
9.466 
10.546 
11.521 

12.426 


8.207 
9.475 
10.556 
11.532 
12.437 


1.81 
2.79 
3.86 
5.04 
6.33 


-0.0313 
-0.0298 
-0.0288 
-0.0280 
-0.0274 


= 1.0 
Zs = 1.5 


5.0 
9.0 
13.0 
17.0 
21.0 


1.30 
1.48 
1.63 
1.77 
1.89 


5.205 
5.933 
6.538 
7.078 
7.576 


5.219 
5.947 
6.552 
7.093 
7.591 


0.74 
1.10 
1.47 
1.87 
2.30 


-0.0396 
-0.0374 
-0.0361 
-0.0351 
-0.0344 



Table 4. The concentration parameters and relative error in the virial mass for a set of possible lensing systems modeled with soft 
truncation versus hard truncation. For soft truncation, the density profile is pi with t = 3c. For hard truncation, the density profile is po 
within the virial radius, and zero outside of the virial radius. The corresponding values of c are both found by integration of the geodesic 
equations. 
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System 


z; 


Zs 


Oe (arc sec) 


RXJ 1347- 1145 


0.45 


0.80 


35 


MS2137+23 


0.313 


1.6 


16 


RCS2319+00A 


0.9 


3.86 


12 



Table 5. Three leasing systems analyzed in Figs. [2] [3] and H System RXJ1347-1145 is from Sahu et al. (1998), system MS2137+23 is 
from Gavazzi et al. (2003) and system RCS2319+00 is from Gilbank et al. (2008). 
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5 10 15 20 

C X 

Figure 1. Density profiles as functions of the scaled distance x = r/rs- The sohd line represents the standard NFW profile x^po{x), 
with a non-integrable fall-off. The other three lines represent the n = 1 truncated profiles x^pi {x) for three different values of the tidal 
parameter r as multiples of the concentration c. All three truncated profiles have integrable fall-ofi' rates. However, the larger the ratio 
r/c, the closer the truncated profile matches the NFW profile within the virial radius, represented by c as a scaled distance. 



40 



• 4 ^ j I I I I ' ' 
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30 - 



-0.01 



T= 3c 



10 - 



-0.02 



^ -0.03 



T= 3c 



. ^ A , 

■ ■ _ * A 4 . . 



-0.04 



0.2 0.4 0.6 
(Mpc) 



0.2 0.4 0.6 
(Mpc) 



Figure 2. Contour in {rs,c) space (left) and relative error between thin lens prediction and actual total mass (right) for the three 

Icnsing systems modeled as truncated distributions pi with r = 3c. Circles correspond to the 35 arc sec arcs in RXJ1347-1145, triangles 
correspond to the 16 arc sec arcs in MS 2137-23, and squares correspond to the 12 arc sec arcs around cluster A of RCS 2319-1-00. 
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40 
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■ * A • • , 



-0.02 



20 - 



T= 2c 



-0.04 



-0.06 



2c 



10 - , . 
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Figure 3. Contour in (j',s,c) sjjacc (left) and relative error between thin lens prediction and actual total mass (right) for three lensing 
systems modeled as as truncated distributions pi with t = 2c. Circles correspond to the 35 arc sec arcs in RXJ1347-1145, triangles 
correspond to the 16 arc sec arcs in MS 2137-23, and squares correspond to the 12 arc sec arcs around cluster A of RCS 2319-1-00. 
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-0.01 



-0.02 



30000000000 



o T= JC 

RXJ1347-1145 



-0.01 



-0.02 



T= JC 

MS2137 + 23 



-0.01 



-0.02 



T= JC 

RCS2319+00A 



-0.03 



10 



15 



-0.03 



10 



15 



-0.03 



10 



15 



Figure 4. Relative error between thin lens prediction and actual total mass for three lensing systems for increasing values of r as 
multiples of the concentration parameter (t = jc). Filled symbols correspond to the truncated model pi and open symbols correspond 
to the non-truncated NFW model po- 



© 0000 RAS, MNRAS 000, 000-000 



